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Multiplicative random processes in (not necessaryly equi- 
librium or steady state) stochastic systems with many degrees 
of freedom lead to Boltzmann distributions when the dynamics 
is expressed in terms of the logarithm of the normalized ele- 
mentary variables. In terms of the original variables this gives 
a power-law distribution. This mechanism implies certain rela- 
tions between the constraints of the system, the power of the 
distribution and the dispersion law of the fluctuations. These 
predictions are validated by Monte Carlo simulations and ex- 
perimental data. We speculate that stochastic multiplicative 
dynamics might be the natural origin for the emergence of crit- 
icality and scale hierarchies without fine-tuning. 

In the last years researchers have found an exceedingly large number of power laws 
in very many natural and artificial (social, economic) systems. 

The emergence of " scaling" properties was considered intriguing as in theoretically 
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known models this is related usually with very special "critical" conditions. In the 
parameter space of typical equilibrium statistical models, critical systems correspond 
to subspaces of measure zero. Yet scaling systems seem to show up in nature much more 
often than this theoretically expected measure zero abundance. This lead researchers 
to coin the term self-organized criticality ( |]], 0). 

In the present note we present a simple yet very general explanation of the emergence 
of power laws. According to our analysis, power-like systems are expected to arise as 
naturally as the Boltzmann distribution. 

In fact we show that for a very large class of systems, their power law distribution is 
in a precise mathematical relation to a Boltzmann distribution when the measurables 
are represented on a logarithmic scale. This analysis implies additional relations which 
are confirmed experimentally. 

Consider a system consisting of a large set of elements % which are characterized each 
by a time-dependent variable Ui{t) (for definiteness one can think of a set of investors 
i = 1, N each owning a wealth Ui or N towns containing each uji people). 

Assume that the typical variations of u are characterized effectively by a multiplica- 
tive stochastic law: 

UJi(t + I) = \uJi(t) (1) 

with v being a stochastic variable with a finite support distribution of probability 7r(A). 

The effective "transition probability" distribution 7r(A) is assumed not to depend 
on i or on the actual value of U{. However, we will see that our conclusions are not 
affected if the shape of 7r(A) varies in time during the process. 

In order to isolate the shape of the distribution of u even for situations in which 
there is an unbounded overall drift of the Wj(t)'s towards infinity, we will work in the 
sequel of the article with the distribution P(w): which fulfills the master equation: 

P(w, t + 1) - P(w, t)= f U(X)P(w/X, t)dX - P(w, t) f II(A)dA (2) 

J A J A 
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where w's are normalized uo's such as to fulfill at each time: 

X>;(*) = / wP(w,t)dw = N (3) 

i.e. in the wealth case one represents actually the relative wealth of each investor. 
Correspondingly, the transition probability distribution n(A) for the new variables is 
related to tt(uj) by a shift in the argument. 

Moreover, one limits from below the allowed values of w > (in the wealth case 
this consists in subsidizing individuals as not to fall below a certain poverty line w ). 
This implies appropriate changes in the transition probability for w^s in the immediate 
neighborhood of w . 

In order to extract the implications of the dynamics (Q) it is convenient to represent 
it on the logarithmic scale in terms of x — Inw and p = ln\. The corresponding 
probability distributions P and II become in the new variables: 

V{x) = e x P(e x ) (4) 

and respectively p(p) = eTI^) . In terms of V,x,p,fi, the master equation (fj) 
becomes: 

V(x,t + 1) — V(x,t) = / p(fx)V(x — /i, t)dfi — V(x, t) / p(p)dp (5) 

J fi J fi 

Not that this equation has the standard form of the master equation for an usual 
Monte Carlo process. 

The iteration of the equation (|) for long time sequences projects upon the eigen- 
mode with the largest eigenvalue of the time evolution operator: 

n p V(x) = J p(n)V{x- n)dn + V{x)(l -J p(p)dp) (6) 

This in turn leads to an asymptotic distribution of V which fulfills an equation of the 
form: 
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/ p{p)V(x - n)dn = AV(x) (7) 
J n 

Ignoring for the moment the boundary and finite size effects, one can easily verify that 
the solution of this equation is: 

V(x) ~ e~ x/T (8) 

with T determined by the condition 

/ eTp(fi)dfi = A (9) 
J [i 

The uniqueness of the solution (|j), @ is insured by the normalization condition (H), 
the positivity of the density distribution V and by the fact that for positive p the left 
hand side in @ is a convex function in ~. A rigorous proof that the equation (0) leads 
to (H|) is given in || and is based on the extremal properties of the G — harmonic 
functions on non-compact groups (in our case the group of translations on R) . 

When one translates back the exponential "Boltzmann" law (Q) in terms of the 
original variables w = e x one gets according ([|) a power-law distribution: 

P(w) ~ w^ 1 - 1 ^ (10) 

If one ignores the departures from (|8p due to the (upper) boundary and finite size 
effects one can use the normalization conditions for the total "wealth", w, eq. (§) 

CI w-^dw = N (11) 



and for the total number of elements: 

poo , 

C / w^-Tdw = N (12) 

J Wo 

in order to express T in terms only of Wq: 

T = 1 — wo (13) 
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This power-law and the above relation [] are excellently confirmed by simulations || in 
various systems for a wide range of wo's and is consistent with experimental data [Q. 

It appears therefore that T is largely independent on the shape of the transition 
probability distribution p(/i) (or n(A)). Physically, an intuitive understanding of this 
result can be achieved by thinking of eq. (|^) in terms of a conservative system in which 
an energy /i can be absorbed or emitted by each degree of freedom i according to the 
("Monte Carlo") emission-absorption probability distribution p(/x). 

The emergence of a Boltzmann distribution is independent on the details of the 
energy exchange mechanism: it is more general than the details of the particular dy- 
namical process leading to it. In fact, even if the process itself is not stationary and 
the "transition probabilities" 11(A) and p(/i) depend on time, the distribution P(w) 
can still converge: modifying during the process (or during a Monte Carlo simulation) 
the interactions from short range to infinite range from 2-body to many-body from 
direct interactions to interactions through the intermediary of a bath or of an "energy 
reservoir" is known || not affect the Boltzmann distribution (|8]). 

One sees therefore that a power law is as natural and robust for a stochastic mul- 
tiplicative process as the Boltzmann law is for an equilibrium statistical mechanics 
system. Far from being an exception and requiring fine tuning or sophisticated self- 
organizing mechanisms, this is the default. 

For our general mechanism to apply to a scaling system, the system has to fulfill the 
effective stochastic multiplicative law (|l]). Yet, the mechanism by which each particular 
system is lead to fulfill ([T]) might differ. For instance in the towns example this might 
be related with interactions between town residents (residents moving upon marrying 



1 For very low values of wq the finite size effects and the upper bound cannot be ignored and 
equation (|l3| ) is modified. The modified relation is confirmed by Monte Carlo simulations too 
|J. In particular for wq = one gets T = oo. 
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somebody from another town, or upon entering a new employee-employer bound). In 
large scale universe structures, 2-body gravitational forces might lead to laws similar 
to ©. 

In a series of papers we have studied in detail theoretically and numerically these 
mechanisms in the context of the stock market ( 0, @, )■ 

Our results turned out to explain in detail the Pareto power law distribution of the 
individual incomes (experimentally documented in 0) as well as the Levy distribution 
in the market prices fluctuations (reported in ]TU|]). 

The Pareto law arises as a consequence of eq. (|T0| ) and of eq. (|J). These equa- 
tions imply that the individual speculative incomes j~i(£) = /(A — l)H(X)wi(t)d\ are 
distributed by a power law flip ) too: 



P(r) = r 1 l -»o (14) 
By a similar argument, one finds that the market price fluctuations induced by individ- 



uals are also distributed according to the fll4|) law. According the generalized central 



limit theorem, a quantity which is a sum of random variables r distributed according to 
a probability distribution r~ a converges to the Levy distribution of characteristic 
exponent 1 — a. 

Our analysis implies therefore (and the experimental available data confirm) that if 
the individual wealth distribution is fitted by a power law of exponent —1 — 1 J w (Fig 
1) then the speculative income distribution is governed by the same law fll4) (Fig 2) 



and the market fluctuations are given by a Levy distribution of characteristic exponent 
-li Fi g (3)- 

These relations are confirmed by the available experimental data (and by the Monte 
Carlo simulation of microscopic representations of the stock market p|) Figs. (1) (2) 
(3) with -a = -j^ = -1.4 JTT|. 



We plan in a future publication to compare with experiment the relations which 
our mechanism predicts between the rate of inflation, the taxation policy and the lower 
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income bound (the poverty line wo). 

In the context of fundamental physics, one may hope that the extension [^] of the 
result (|8|)- ([T0|) to G — harmonic distributions on general non-compact (Weyl gauge 
groups G might lead microscopic models naturally, without fine tuning, to criticality 
and scale hierarchies. 

Such discrete (lattice gauge) theories with non-compact group might provide a uni- 
fied context for treating renormalization theory and time: the continuous re-scaling (|3|) 
of the "running to infinity" degrees of freedom u>j ([TJ) suggests A* as the microscopic 



stochastic origin of both time flow and renormalization flow (with the "extremal" [12 
distribution P as the fixed point). 
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Figure Captions 
Figure 1: 

The distribution of wealth (log-log scale). The solid line represents actual data 
from Great Britain (Source: Inland Revenue Statistics, 1970) Dots represent the wealth 
distribution in our simulations. The dashed line is a fit by a power law distribution 
with a slope of - 2.40 (a = 1.4). 

Figure 2: 

Empirical distribution of income (log- log scale). Data from Great Britain (Source: 
National Income and Expenditure 1970). Dashed line represents fit by a power law 
distribution with slope of - 2.34. 

Figure 3: 

Distribution of returns on the stock (semi-logarithmic scale). The solid line rep- 
resents the Levy distribution with exponent a = 1.40 ( scale factor 0.00375). Dots 
represent distribution in simulation. Diamonds represent the empirical return distri- 
bution for the S&P 500 index during 1984 - 1989 as reported by Mantegna Sz Stanley 
fT0|j . The dashed line represents the Gaussian distribution with the empirical standard 
deviation (a = 0.05) . 
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